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We demonstrate an in-line digital holographic microscopy using a consumer scanner. The consumer 
scanner can scan an image with 4,800 dpi. The pixel pitch is approximately 5.29 jum. The system using a 
consumer scanner has a simple structure, compared with synthetic aperture digital holography using a 
camera mounted on a two-dimensional moving stage. In this demonstration, we captured an in-line 
hologram with 23, 602 X 18, 023 pixels (—0.43 gigapixels). The physical size of the scanned hologram is 
approximately 124 mm X 95 mm. In addition, to accelerate the reconstruction time of the gigapixel 
hologram and decrease the amount of memory for the reconstruction, we applied the band-limited 
double-step Fresnel diffraction to the reconstruction. 



Digital holographic microscopy (DHM) captures a hologram with an electronic device such as a CMOS or 
CCD camera, and the captured hologram is reconstructed on a computer using diffraction calculation 1 ' 2 . 
In order to increase the field- of- view, lateral and depth resolving powers of the reconstructed image, we 
need to capture a large hologram, for example, the amount of gigapixels achieved in recent researches 3 " 6 . There are 
several methods available to capture a gigapixel hologram. In astronomy, a single CCD device with over 0.1 
gigapixels has been achieved; however such a large area CCD is expensive. Another method for acquiring a 
gigapixel hologram is synthetic aperture digital holography using a camera mounted on a two-dimensional 
moving stage (or moving a reference light). Reference 5 reported a gigapixel digital holography using a typical 
color CCD that has Bayer color filter, and multiple holograms obtained by changing illumination angles of a 
reference light. A hologram captured with such a color CCD sensor and green light source lacks the pixels 
corresponding to the red and blue due to the Bayer filter, that is the effective pixel number is decreased. In order 
to increase the pixel number, Ref. 5 proposed ideal interpolation technique of the lacked pixels by using adjacent 
pixels. Recently, gigapixel microscopy using a consumer scanner has been proposed 7 . The approach is excellent 
because the microscopy has a wide field-of-view, low-cost and simple structure; however, it cannot observe a 
sample in the depth direction. 

In this paper, we demonstrate an in-line DHM using a consumer scanner, inherently observing a sample in 
three-dimensions. In-line DHM 8 ' 9 is capable of obtaining a hologram without using beam splitters and mirrors. 
The consumer scanner can maximally scan an image with 4,800 dpi. The pixel pitch is approximately 5.29 fim. 
The DHM system using a consumer scanner has a simple structure, compared with synthetic aperture digital 
holography. In this demonstration, we captured an in-line hologram with 23, 602 X 18, 023 pixels (—0.43 
gigapixels). The physical size of the scanned hologram is approximately 124 mm X 95 mm. In addition, to 
accelerate the reconstruction time of the gigapixel hologram and decrease the memory usage for the reconstruc- 
tion, we applied the band-limited double-step Fresnel diffraction (BL-DSF) 10 to the reconstruction. 

Results 

We show a hologram of the USAF1951 test target captured by the DHM and the reconstructed image. Figure 1 
shows the hologram with 23, 602 X 18, 023 pixels (0.43 gigapixels) and the hologram is sampled at about 5.29 fim. 
The hologram is captured under the condition that the distance between the scanner and the sample is 50 cm and 
the distance between the sample and the light source is 15 cm. Because the number of pixels in the hologram is 
much larger than the resolution of a display (1,920 X 1,080), the decimated hologram by image processing is 
displayed. The inset shows the raw hologram in a part of the hologram (in the red square). We can observe the 
interference fringe of the hologram. Although the scanner can obtain a hologram maximally in 16bit/pixel, we 
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Figure 1 | Hologram of USAF1951 test target with 23, 602 X 18, 023 pixels and the hologram is sampled at about 5.29 /an. 



were unable to recognize the difference between the reconstructed 
images in 8bit and 16bit/pixel. Therefore, we used a hologram cap- 
tured by 8bit/pixel in terms of the calculation time and the memory 
usage. 

The reconstructed images are shown in Figure 2. The obser- 
vational area of Fig. 2(a) is about 22 X 29 mm 2 . Figure 2(b) shows 
the details of the red square in Fig. 2(a). Figure 2(c) shows the red 
square in Fig. 2(b) in greater detail. We can observe a reconstructed 
image in the resolving power of about 8.8 /urn. The movie shows the 
observation of the reconstructed image as magnifying the image. 

In a large hologram, the calculation time and memory usage for 
the reconstruction are important issues. We compared the calcula- 
tion time and memory usage in the reconstruction of the hologram 
(Fig. 1) using the angular spectrum method (ASM) 2 and BL-DSF 10 . In 
the calculation, we used Linux (64-bit) as the operating system and 
Intel Core i7-2600S as the CPU, the PC has the memory amount of 16 
Gbytes. In ASM, the calculation time and the memory usage were 
about 355 seconds and 12.6 Gbytes, while the calculation time and 
the memory usage of BL-DSF were about 177 seconds and 3.2 
Gbytes, respectively. The peak signal-to-noise ratio (PSNR) between 



the reconstructed images of ASM and BL-DSF is over 30 dB. ASM is 
widely used in digital holography community as the reconstruction 
method. Therefore, we use reconstruction results obtained by ASM 
as the criteria to other reconstruction method, namely BL-DSF. 
Generally speaking, the difference of two images is very few when 
PSNR is over 30 dB. 

Figure 3 shows a hologram that, records an ant and water-flea 
placed at 30 cm and 50 cm from the scanner, respectively, and the 
reconstructed images. The number of pixels of the hologram is 23, 
602 X 18, 023. When focusing on the water-flea, the ant is unfocused. 
While, when focusing on the ant, the water-flea is unfocused. 

Discussion 

Figure 4 shows the outline and photograph of the in-line DHM using 
a consumer scanner, respectively. In-line DHM 8,9 is capable of 
obtaining a hologram without using beam splitters and mirrors. 
Samples are placed between a light source and the hologram. The 
diffracted and undiffracted lights by the samples are regarded as 
the object light and the reference light, respectively, then, the 




(c) 



Figure 2 | Reconstructed images from the hologram of USAF1951 with 23, 602 X 18, 023 pixels. (See supplementary video). 
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Figure 3 | Hologram that, records an ant and water-flea placed at 30 cm and 50 cm from the scanner, and the reconstructed images. 



interference fringe of these lights generates a hologram. In-line DHM 
can control the area of the hologram and the magnification of the 
object light just by the location of the light source or the numerical 
aperture (NA) of an objective lens. If using a holographic setup that 
requires beam splitters, lens and mirrors, a large aperture is needed 
because the scanner has a large aperture; therefore, the in-line DHM 
is suitable for gigapixel holograms. 

We used a fiber-output laser with a wavelength of 405 nm. 
Although the output of the laser has a spherical wave, we used an 
objective lens with a magnification of X10 and an NA of 0.25, to 
expand the angle of spread of the fiber- output laser. We placed 
samples between the objective lens and the scanner. The scanner 
captures an in-line hologram by moving the image sensor of the 
scanner. 

As shown in Fig. 5, consumer scanners are mainly categorized by 
two types: "CCD" and "Contact Image Sensors (CIS)" scanners. 
Regarding CCD scanners (Fig. 5(a)), they have a two-dimensional 
CCD sensor whose size is smaller than the scan surface (cover glass); 
therefore, a reduction optical system composed of some mirrors and 
a lens is required to reduce the image on the scan surface to the 
CCD. And CCD sensors generally have a color filer for scanning 
color images; however, the color filter disturbs the capturing of 
holograms 5 . 

In contrast, CIS scanners (Fig. 5(b)) have a simpler structure than 
CCD scanners because CIS sensors are one- dimensional devices 
whose size is the same as the scan surface; therefore, reduction of 



the optical system in CCD scanners is not required. In addition, 
because CIS sensors are one dimensional, the electronic circuit is 
simpler than CCD sensors; therefore, the CIS scanners can capture 
an image with 16bit/pixel maximally, while most CCD scanners 
capture an image only with 8bit or 12bit/pixel maximally. CIS sen- 
sors also do not need a color filter unlike CCS scanner because CIS 
sensors capture a color image by switching RGB light sources in 
time-divisions; therefore, CIS scanners are not adversely affected 
by the color filter. Thus, we adopted a CIS scanner in this research 
because CIS scanners are suitable for hologram recording. 

We used "CanoScan LiDE 210" made by Canon as the scanner. 
The scanner has a CIS image sensor and can capture an image with 
the same resolutions (4,800 dpi) in the horizontal and vertical direc- 
tions. We used "ScanGear (32 -bit version)" as the scanner software. 
The GRIN lens attached to the CIS sensor is used for the imaging of 
scanning plane. However, the GRIN lens hamper a hologram record- 
ing. Therefore, we removed the GRIN lens. We covered the light 
source by a black tape because of cutting off the light. The scanner 
can maximally scan an A4 size image (297 mm X 210 mm) with 
4,800 dpi, theoretically achieving a resolution of 56, 144 X 39, 698 ~ 
2.22 gigapixels. However, we could capture a hologram with 23, 602 
X 18, 023 pixels because we might use 32-bit version of the software. 

Methods 

The reconstruction of gigapixel holograms is very time-consuming and requires 
huge memory usage. We adopted BL-DSF 10 , which is an effective way to obtain a 
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Figure 4 | In-line digital holographic microscopy using a consumer scanner. (Left) Outline of the system. (Right) Photograph of the system. 
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reconstructed image from gigapixel holograms. In Fourier optics, diffraction calcu- 
lations are categorized by two forms: the first is convolution-based diffraction and the 
second is Fourier transform-based diffraction. Here, we show ASM as an example of 
convolution-based diffraction. ASM is expressed as follows: 



u 2 (x 2 ,y 2 ) =FFT 



FFT[ Wl ( Xl , yi )] exp I - InizJ 1 / X 2 -ft -ffj 



(1) 



where X is the wavelength of light, the operators FFT[*] and FFT _1 [*] are the fast 
Fourier and inverse fast Fourier transform respectively, u x (x ly yi) and u 2 (x 2 , y 2 ) 
indicate a source and destination planes, (f x ,f y ) is the coordinate on the frequency 
domain and z is the propagation distance. The merit of convolution-based diffraction 
is that the sampling rates on the source and destination planes are the same; however, 
the demerit is the need to expand the source and destination planes by zero-padding 
to avoid noise. The form of ASM is the convolution using FFT. Generally speaking, 
the numerical calculation of the convolution using FFT that has circular property 
requires zero-padding to be the size of 2N X 2N, where N is the horizontal and vertical 
pixel numbers of the hologram, during the calculation. If we do not use the zero- 
padding, the reconstructed image will incur noise by wraparound due to the circular 
property. The zero-padding can avoid the wraparound. The usage memory and 
calculation cost of ASM are proportional to 4N 2 and 4A/ 2 log 4AT, respectively. It takes 
large memory usage and long calculation time. 

To overcome this problem, double-step Fresnel diffraction (DSF) has been pro- 
posed 11 . It calculates the light propagation between the source to destination planes by 
twice the Fourier transform-based diffraction, via a virtual plane (x v ,y v ). DSF does not 
need zero-padding because DSF is based on Fourier transform-based diffraction. In 
addition, although most Fourier transform-based diffractions change the sampling 
rates on the source and destination plane, DSF does not change them. Due to the 
original DSF incurring the aliasing noise under certain conditions, we proposed BL- 
DSF introducing the rectangular function for band-limitation to the original DSF. BL- 
DSF is expressed as follows: 



2 (m 2 ,n 2 ) = C Z2 VVV^ 



exp \ v /vJ Rect 
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(2) 



where z x is the propagation distance from the source plane to the virtual plane, z 2 is 
the propagation distance from the virtual plane to the destination plane, 



C Zl = exp ^- — (x 2 2 +yl) j • The operator FFT^" (z) means the forward FFT when the 

sign of z is positive and the inverse type when it is negative. See Ref. 10 how to 
determine the band-limiting area {x™ 11 *, y™"*). The form of BL-DSF is not the con- 
volution. The usage memory and calculation cost of BL-DSF are only proportional to 
N 2 and N 2 log N, respectively. BL-DSF was implemented in our open-source wave 
optics library, CWO+ + 12 . 
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